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"ti ' Abstract 

We consider eigenvalue problems in quantum mechanics in one dimension. Hamil- 
Q" 1 tonians contain a class of double well potential terms, x 6 + ax 2 , for example . The 

space coordinate is continued to a complex plane and the connection problem of 
fundamental system of solutions is considered. A hidden U q (gl(2\l)) structure arises 
in "fusion relations" of Stokes multipliers. With this observation, we derive coupled 
nonlinear integral equations which characterize the spectral properties of both ±a 
potentials simultaneously. 
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1 Introduction 



The eigenvalue problem of a one-body ID Schrodinger operator is the most fundamental 
subject in quantum mechanics. Still, it provides vivid materials of research. 

Besides a few exceptions where eigenvalues and wavefunctions are obtainable explicitly, 
one may employ several tools for analysis, e.g., the perturbation theory, the variational 
approach and so on. Among them, the exact WKB method |1[]-[BJ is unique in the sense 
that it provides non-perturbative information on the analytical structure of wavefunctions 
and spectral properties. We analytically continue x, original coordinate variable, to a 
complex number. The whole complex plane is divided into several sectors. In each sectors, 
there are two linearly independent solutions , as Schrodinger operator is the 2nd order 
differential operator. They are referred to as the fundamental system of solutions (FSS) 
in the sector fl8|| . The relations among FSS in different sectors are central issue in the 



connection problem. The importance of the problem and consequently Stokes multipliers, 
in the WKB problem has been deeply recognized and emphasized in early '80, especially 
in !• 

Recently, a remarkable link has been established among the spectral determinants of 
a ID Schrodinger operator associated with the anharmonic oscillator, transfer matrices 



and Q operators in CFT possessing U q (sl(2)) ||-[]12|. Here the spectral determinants 
imply D(E) = n% eigenvalue (1 — if") an d hs generalizations. A curious interplay between 
D(E) and generalized Stokes multipliers is also found |], [10D . In view of solvable models, 



a striking fact is that they share the same functional relations with transfer matrices in 
the fusion hierarchy possessing U q (sl{2)) |I| [13], [14], p~5]| . This allows for applications of 



the strong machineries in the study of solvable models [13]- [27] to the studies of Stokes 
multipliers, spectral determinants and so on. Several results have been explicitly obtained 
for the anharmonic oscillator problem, and are extended to higher differential analogues 



1,0 



In this note, we consider an anharmonic oscillator perturbed by a lower power poten- 
tial term. It belongs to a class of potentials discussed generally in || [7J with the exact 
resolution method. To be precise, we consider the eigenvalue problem, 

H(x,a)* k (x) = (-^+^ + a x M4 )^(x) = EA(4 (1) 

Throughout this report we set h = 1 and M > 1. 

The spectral problem concerning this Hamiltonian turns out to be in a category to 
which one can apply the tools in solvable models. 

The sign of a seems to be crucial if one considers the operator ([[]) on the real axis. 



FIG.l 



We will not expect much difference from the "pure" anharmonic oscillator when a > 0, 
while we do expect change for a < as the potential develops the double well. 
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It will be shown, however, that the negative a and the positive a problems are not 
separable when we discuss the global connection problem. Roughly speaking, the negative 
a problem is coupled to the positive a problem by crossing a border line of neighboring 
sectors and vice versa. See §2 for precise arguments. It may be then reasonable to consider 
a two-fold connection problem (crossing two adjacent lines) , or more generally, relations 
between sectors separated by even multiples of border lines. Some of the Stokes multipliers, 
in the generalized connection problem, possess expressions corresponding to the eigenval- 



ues of the (fusion) transfer matrices of the 3 state Perk-Schulz (PS) model ^9|, |30| of 
which underlying symmetry is U q (gl(2\l)). Others can not be directly equated with the 
(fusion) transfer matrices but have relations with the the 3 state PS model as well. Thus 
we conclude that the perturbation ax M ~ x breaks the U g (sl(2)) symmetry of the "pure" 
anharmonic oscillator but it brings the new symmetry U q (gl(2\l)). The deformation pa- 
rameter q is related to the exponent of the perturbation by q = exp(ij^). Through these 
findings, we can derive the nonlinear integral equations (NLIE) which characterize the en- 
ergy levels of both the negative a problem and the positive a problem simultaneously. 

The paper is organized as follows. In the next section, we will explore symmetries 
of solutions to (H). The precise definition of sectors is given. The connection problem is 
addressed in section 3. Certain components in fusion Stokes matrices are identified with 
eigenvalues of fusion transfer matrices associated to U q (gl(2\l)) . The spectral determinant 
is explicitly parameterized by FSS in a sector. The coupled NLIE are then derived in sec- 
tion^ which determine energy levels. We will perform analytical and numerical checks on 
the consistency of our result in section 5. Section 6 is devoted to summary and discussions 
on open problems. 

2 Asymptotic Expansion and Symmetry of solutions 

Let 4>(x, a, E) be an entire function of (x, a, E) and a solution to TC(x, a)4>(x, a, E) = 
E<p(x,a,E). 

The solution, which decays exponentially at x — > oo, is of primary interest. By em- 



ploying the argument in |31j, we immediately find its asymptotic behavior, 



t m+i 

<P(x,a,E) ~ x-^expl-^), (2) 

r M+l 

d x <f>(x,a,E) ~ X "/W2e^(___). ( 3 ) 

The validity of the above expansion is not restricted to the real axis, but extends to the 
wedge, |argx| < 2 ^ +2 ]T3"| [JT|. 



The second order linear differential equation admits another independent solution. To 
specify it, or to deal with the global problem, it is convenient to extend x to the complex 
plane as mentioned in introduction. Then, as in the case of a = 0, the solution exhibits a 
symmetry by rotating the complex x plane by a specific angle. 
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The direct calculation proves the following. 



Theorem 1 Let <p(x, a, E) be the above solution and q 



Af+l- 



Then 



4>(q x, q M+ a, q E) is also the solution to the differential equation, H(x, a)(j) = Ecj). 



This is the desired second solution which grows exponentially on the positive real axis: 



x 



-M/2+a/2 



exp( 



M+l 



for x — ► oo We note that q 



M+l 



a 



—a. This deserves an attention. 



As mentioned in introduction, the potential assumes the completely different structure for 
a positive and a negative on the real axis. The rotation in the complex x plane by angle 
however, couples these two problems. Thus we shall treat the Hamiltonians with 



M+l ■ 



±a simultaneously. Similar pairing of differential equations is found for the positive and 



the negative angular momentum terms in a class of 3rd order differential equations [16 



it 



This observation is crucial in our approach and can be generalized further. To state 
we prepare some notations. 
Hereafter a always takes a non- negative real value. By TC^ e \x, a), we mean the 
Schrodinger operator, 

d 2 



dx 2 



+ x 2M + eax M ~ 



where e = ±1. 

Let Sk be a sector in the plane satisfying 



|argx 



< 



71 



M+V ~ 2M + 2 



FIG.2 



The FSS depends on the sector. We conveniently define 

yf ■= 7= (j){xq 3 ,ea,q 3 E). 



Theorem 2 For the 1-L^ e \x,a), the FSS in the sector Sj is given by (yf 3 ^, ) where 



For a = case, this has been argued in |13|, [10| . It is easily checked that is the 
sub-dominant solution in Sj] it tends to zero as x tends to infinity along in any direction 
in the sector. 

In the next section, we consider the global connection problem of these FSS in the 
complex x plane. 
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3 Fusion Stokes multipliers, U q (gl(2\l)) structure and 
Spectral Determinants 

We introduce the Wronskian matrix 

\d x yy, d x y) +1 'J 

and the Wronskian Wff' := det$^(x). 

The linear dependence of the solution can be easily verified by evaluating the Wron- 
skian at Sj+1/2 using the asymptotic expansions (0) and (0) . The present normalization 
yields W fc (e) = 1. 

Let A4jl be the Stokes matrix connecting the Wronskian matrices (a:) and 3>j-+i (x), 

*f\x) = *$(x)M$. (5) 
It permits an explicit parameterization 

"53 - g; J) • (•) 

where r- is referred to as the Stokes multiplier. We have two remarks. First, the (1,1) 
element is the function of a and q^E. We omit the dependency on a. The dependency 
on E is indicated by the index j. Second, the (2,1) element,— 1, is a consequence of the 
present normalization of the Wronskian. 

To be more specific, we consider the operator TiS + '(x,ot) and start from the positive 
real axis (or more generally S ). The initial FSS is (y^\y[ •*). 



FIG.3 



A linear relation follows from 
neighboring sector Si, 



between this FSS and (y[ , y^ ), the FSS at the 



r (+y-) 

r ill 



(7) 



Similarly, the FSS in 5 2 is linked to the FSS in Si by 



v{- ] 



r (-y+) 

'0 4/2 



Judging from the upper indices which indicate the corresponding signs of a, it may 
be natural to introduce a generalized Stokes matrix A4q2 connecting FSS (yo + \ y[~^) and 



(y 



<+) 



. It is simply obtained by the matrix multiplication, 
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- /Vl 0,2 — Jvl l,l - /vt 0,l 



r (-) r (+) _ j r (-r 



—T, 



(+) 
) 



-1 



(9) 



Equations (|7|) and (H) yield r's in terms of y's. The (1,1) component in (||), hereafter 
denoted by Ti tl (E), is then represented in terms of y as, 



Ti.i^H^-W-i 



yi> +) , + yT } 



y^ 



(10) 



The dependence of in the rhs is implicitly indicated by indices of y. We will comment 
on this representation in terms of a solvable model later. 

There is another expression using both y's and dy's. This form is of practical use in the 
following generalization. By applying the Cramer formula to (|5j), we immediately obtain 



(<0 



d x yf, d x yfl 2 



Note that we use the fact that the Wronskian is normalized to be unity. The (1,1) entries 
in (|9|) is then given by 

d x yo + \ d x y ( 3 } 

One can further generalize the above result. Naturally, the "fusion" Stokes matrix 
M.^2k ls defined which relates FSS of Sj to Sj + 2k- Explicitly, it is given by 



M 



(+) 

0,2fc 



V 



„(+) Jr) 

yo > y2k+i 
d x y , d x y 2k+1 

<9xyo j dxVik 



y(-) J') 

i/1 5 i/2fe+l 

<9x2/i j d x y 2k+1 
d x y[ \ d x y { 2 l ] 



(12) 



for j = 0. We can prove the above using the induction on k most easily. Similar formula 
holds for _Mq 2k by replacing all upper indices + <-> — . 

We are now ready to relate an entry in a fusion Stokes matrix to the spectral deter- 
minant. Hereafter we assume M = 2m — 1. It follows from the above argument that 



2m) 



(13) 



FIG.4 
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y Q ( y-zm ) stands for the subdominant solution on the positive (negative) real axis. 
They tend to zero asymptotically in their proper region, being appropriate basis for the 
eigenfunction. The eq.(|T3D tells, however, that y^ is combined to both y$ and y[ ■* by 
rotating the complex plane by — ir, 



vt2 = {c iy ^ ] + c 2 y i f ) )/detM^ m 



This is an obstacle in constructing an eigenfunction defined on the whole real axis. The 
prescription is to demand that the coefficient of y[ ^ must vanish if E is an eigenvalue. 
Consequently, it is proportional to the spectral determinant. 

The coefficient is essentially equal to the (2,1) component of M.Q2 m i anc ^ ^ reads in 
terms of the original function as 

g Q(m+1) <j>{x,a,E), <j>(-x,a,E) 
2 d x (f>(x,a,E), d x (f>(-x,a,E) 

The x dependencies are spurious as the entities are products of Stokes multipliers 
which are obviously x independent. We adopt the simplest choice x — 0. The coefficient is 
now proportional to 0(0, a, E)d x (f>(x, a, E)\ x= q. Thus we conclude that for an eigenvalue 

0(O,a,£j +) ) = or d x <f>{x, a, E$ +) )\ x=0 = (14) 



E\ +) of H(+)(i,a) 



must hold. 

We can repeat the same argument starting from TiS~\x, a) on the positive real axis. 
The above observation may lead to the identification 

<f>(0,ea,E) ~ DW(E):=IJ(1-E«), (15) 

i 

d x( p(x,ea,E)\ x=0 ~ d¥(E):=1[(1-EV). (16) 

3 

The lower sign signifies the parity: the positive parity means a contribution from 
symmetric wave function. The product must be taken over eigenvalues with the corre- 
sponding parity. The total set eigenvalues {E^} of H^(x,a) consists of two subsets, 

{Ef} ={EQ} U {E®j} and D^(E) = D^(E)D { 1\e). 

We comment on the relation of the present result to an existing solved model. T]_^(E) 
in (|10D can be represented, utilizing ([16|) , as 



D { -\E) 2a D M (E)Di-\q«E) +1 Dt\g e E) 
D { +\q*E) D?\q*E)Dt\q 2 E) D ( l\q^E) 



where we safely choose x = in the rhs. 
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The above expression has similarity to the dressed vacuum form (DVF) of the (un- 
fused) transfer matrix for the 3-state PS model with grading (+, — , +). The latter can be 



found in eq.(3.1) and (3.2) of []30], The spectral parameter v corresponds to energy in 
the Schrodinger operator, precisely, E = exp(-^y). 

The spectral determinants have the following identification to the eigenvalues of Bax- 
ter's Q operators, 

D^\E) = Q 2 ( v + J -^i) 



Those with positive parity may be identified with the second solutions of Baxter's Q 
operators. 

The scalar factors (vacuum expectation values) f a (x), g a (x) in |32| depend on the choice 



of the quantum space. We assume that the present quantum space space gives the simple 
scalars as in (|17|). In this sense, Tx,i(E) exhibits the hidden U q (gl(2\l)) symmetry behind 
the present Schrodinger operator, just as in the U q (sl(2)) symmetry for a = problem. 
This coincidence can be observed further. We have checked up to certain value of k that 
the (1,1) element and the (2,2) element of M-Q^k coincide with DVF of symmetric fusion 

transfer matrices A^ and — Ajjjj^ in [^], respectively. The interpretation of the (1,2) and 



the (2,1) element, in terms of fusion transfer matrices, is still an open problem. 
One can adopt another description of T\ )1 . The (2,1) component of (|5|) results 

Proceeding as above, we arrive at, 



T (F)-n^ I n 2* D { + + \E)D { + -\Q 6 E) . D^^E) 

D { +\q±E) D { +\qiE)D { -\q*E) D { -\q 2 E) 

In the next section, we determine the energy levels by utilizing the above results. 
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4 Nonlinear Integral equations for eigenvalue prob- 
lem 



The Bethe ansatz equations follow from the pole-free property of T\^(E) on the real E 
axis, 



(20) 



7 i^V 2 <j) " ^ +) (g- 2 <j^ ^ 



To their analysis, we apply the strong machinery in solvable models, the method of nonlin- 
ear integral equations. Hereafter we shall confine ourselves to the case < a < M where 
energies are nonnegative. The simple pattern of energy spectrum permits the following 
simple-mind choice of auxiliary functions, 



4\E) := r- ^f" q :> (22) 



D^\ q -*E) 

Af(E) := 1 + 4\E). (23) 



Thus 



Af(Efj)=Q. (24) 



Remember that e(= ±1) represents the signature of the perturbation while e'(= ±1) 
denotes the parity. 

In addition, we need some inputs about the asymptotic behaviors from the WKB 
method. Fortunately, they are already available as the existence of lower power term 
does not alter them. 



\nD [ l ] 



(E) 



b E 



00 
i 

/' 

b 



a 



-Ey 



2vr(j + -), 
M+l 



a 



\E\ — > oo, |arg(— E)\ < n 

— > oo 
bo 



(25) 
(26) 



2M 
MT(- + — ) 



2 sin fin ' 



There might be several routes to reach nonlinear integral equations among 
A$y Here we choose the quickest way || [1 



21], E3] which fully exploits the fact that 



and 
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zeros of D±\E) are on the positive real 9 axis. In addition, we assume that there are no 
zeros of A±\E) inside the narrow strip including the positive real axis other than those 
from zeros of D±\e). 
Apparently we have 

fei/^ (eat — e')-n v-^ ^, E . 

3 e',3 

The above assumption allows the representation of the summation part by an integral 
over contour Ce which surrounds the positive real axis counterclockwise, 



log4 e) (£) = ^ ^ + ±f dE>F(^)d E , logA^iE'). (27) 



(cn-c')n. 1 r E 

>c E 

For convenience, we introduce a variable 9 || by 



E = exp(#//i)/V 2 u= (2M + 2)~^/r(— ) 



which originates from the matching condition of the WKB result (|25|) and the Q— operator 
analysis |23|, |33| . 

Let o[, (9), (0) be auxiliary functions defined in (P2p, ( |2"B| ) regarded as functions of 
9. Then eq.([27]) reads, 

loga?(0) = {Ea M ^ i + ^ jf d9'G(9 - 9')d 6 , log^V) 

G{0) = log^ ^f 1 + 

smh (MTT " z mTi) 

where Ce encircles the whole real axis counterclockwise. 

For the reason which will be supplemented, we shall keep (9) in the lower half plane 
but use 1 / a*? (9) in the upper half plane. 

This requirement modifies the above expression as 

loga?(0) = fc^^i * f°° de G(9-9' + iO) i og a { e 7 e) (9'- iO) 

M + 1 27U y,^ 

/ d9'd e G{9 -9' + iO) \og% { J e) {9' - iO) 

where 6 1 is assumed to possess small negative imaginary part. The property (a^\9)) 
-7-r — is employed in the above transformation. 
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We solve (p§[) in terms of loga^(#) to reach the final expression of NLIE, 



/oo 
#i(0 - 9' + zO) In 21^(0' - i0)d6' 
-oo 

/oo 
K 2 (9 -9' + iO) In 2^(0' - »O)d0'}. 
-co 



(29) 



The kernel functions read 

K^9) = - 



1 

2^ 



„iw6 blml 2M 



Few remarks are in order. 



— oo 

J /"OO 

27T j _ N 



sinh ttw sinh ^ 



Awd k 



smh 2A / sinh 



2Af 



sinh 7rw sinh ^ 



1. As a consequence of connection rules, the integral equations are coupled for auxiliary 
functions related to the positive and the negative coefficient of x M ~ x . 

2. On the other hand, equations with different parities are decoupled. 

3. The constants are determined from the consistency by putting 9 — ► — oo. The a 
dependence is only summarized in these constants. 

4. The first term in the rhs is determined so that we recover the result from the WKB 
method Q25D by dropping contributions of integrals. Clearly, (9) is bounded in the 
upper-half plane. This explains our choice of appropriate half planes for auxiliary 
functions. 

The eigenvalues {-E'j > ±} are evaluated by 

In a ( ±\9\%) = {2j + and E$% = exp(^//i)/z/ 2 . 

More explicitly 

i&oi/-V# = (2j + 1 - e'~ + e— )tt 
2 ° v J 2 2M J 

/OO POO 
- 0' + iO) lnSt^tf')^' + / #2(0$ - 0' + iO) ln2l<r e) (0')^'} 
-oo •/ — oo 

(J > 0). (30) 

We present examples of numerical solutions to in Fig. 5. The real and the imagi- 
nary parts of In 21+ are depicted for M — 3, a — 1. 



FIG.5 
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5 benchmarks 



We shall check the nonlinear integral equations analytically for limiting cases and numer- 
ically. 

(1) a = case 

By putting, (E) = a^T (E) the coupled NLIE reduce to an identical integral equation. 
Immediately seen, the result coincides with the nonlinear integral equation in ||. 

(2) a = M case 

In this case, we have a duality in energy spectra; {Ej + ^} coincide with {Ej} , except 
for Eq~^ = in the latter. This degeneracy can be easily explained by the following 



representation of the Hamiltonians |34|, p5 



H { -\x,a = M) = Dtp 
H {+) (x,a = M) = VV ] 

_ 1 d 
T> = - — ix 

i ax 



M 



Once an eigenvector Tl^~\x, M)ipj ' = Ej ^ is found, we can construct the eigenvector 

for 7i^ + \x,a) with the same energy by V^-i := T*^] ■ Only the exception is the j — 

case where T>ipjJ = 0. It is interesting that the asymptotic form (0) from the WKB type 
argument is exact for all x in this case. [] 

The above facts can be also verified from (p9|). Note that the rhs can be treated as 
the mod 2ni quantity. Then the choice a = M leads to the same coupled equations under 
identifications *-> ~\0), a M (6) <-> \0). This explains the de generacy of the 

spectra as it consists both from the negative and the positive parity contributions. The zero 
energy case must be treated more separately. By choosing j — 0, e = — e' — — l,a — M, 
we find the first term in lhs of (p0|) is null. So the first order approximation is 9±q — — oo. 
Actually this is exact as we determine the constant terms so that NLIE are consistent in 
6 — ► — oo. See remark 3 after ( p9|) . This solution gives the missing energy 0. 

Finally we present the preliminary numerical results for M = 3. 

Table 1 shows the results from the IMSL package (dsleig.f), the (naive) WKB method 
and those obtained by solving the nonlinear equations. The agreement is not yet precise 
enough (typically 3-4 digits). Some implement is still in need for the numerical accuracy. 
Nevertheless, the NLIE data already show much improvement from the (naive) WKB 
results. 

By the (naive) WKB method, we mean a self-consistent determination of E^ by 

Xo 



\p\dx = I \J E^f 1 — x 6 — eax 2 dx = (j + —)it (31) 



1 I thank V.V. Bazhanov, R.J. Baxter and B. Nienhuis for pointing out the explicit eigenfunction for 
i/jq ' right after my talk. 
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where Ej = 0. Particularly, for the value with asterisk, this method has 

subtlety. Immediately seen, Eq = 0, xq = 2 1//4 is a formal solution to ([H]) for j — 0, e = 
— ,a — 2. It however involves an isolated turning point of the 2nd order at the origin if 
E^ = 0, which spoils the simple application of the condition (|31~|). The value with o has 
similar difficultly We however skip further discussion on the validity on the (naive) WKB 
method as it is out of the present subject. 

Summarizing, we check the consistency of ( p9|) in some limiting cases and by numerical 
methods. 

6 Summary and Discussion 

In this report, the eigenvalue problem has been addressed for the ID quantum sys- 
tems of which Hamiltonians include double well potentials. We have successfully derived 
the coupled NLIE which determine energy levels of the systems with potential terms of 
±ax M ~ ! + x 2M at the same time. 

The essence of our strategy is to utilize the following correspondences between ID 
quantum mechanics and 1+1 D solvable models, 

energy spectral parameter 

Stokes multipliers -<=/- transfer matrices 
eigenfunctions or derivatives at x = -<=/- vacuum expectation values of Q operators 

We are then entitled to apply the strong machinery of the latter developed since 
Baxter's revolution. 

There are several open questions. 
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1. In this report we confine ourselves to the simplest case a < M. For a > M, the ex- 
istence of negative eigenvalues ruins the analyticity assumptions on auxiliary func- 
tions. Still, formal expressions of NLIE are possible which are similar to excited 
states TBA equations. The integration contour is, however, not so simple as de- 
scribed here. This is an apparent drawback in actual numerical investigations. The 
clever choice of auxiliary functions may be desired. 

2. The understanding is lacking on the intrinsic reason why afflne symmetry like 
U q (sl(2)) or U q (gl(2\l)) comes into play in this simple ID quantum mechanical 
model. 

3. This is somewhat related to the above, but is the most intriguing question. Where 
is the Yang-Baxter equation in the ID Schrodinger operator problem?. 

Once this is known, the fusion hierarchy, useful in the present study, is a mere 
corollary of it. 

We hope to answer these in the future publication. 
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FIGURE CAPTIONS 



Fig 1. An anharmonic oscillator perturbed by a positive or a negative perturbation term. 
Fig 2. The complex plane is divided into sectors. So and Si are indicated as examples. 
Fig 3. The FSS in S and S\ are related by the matrix M^}. 

Fig 4. The connection of FSS on the negative and the positive real axis is accomplished 
by M { %. 

Fig 5. Left: the real part of In 21^, Right: the imaginary part of In 21+^. 

TABLE CAPTIONS 

Table 1. First two energy levels calculated by the IMSL library, the (naive) WKB method 
and the NLIE method. We choose M = 3 and adopt various a. See asterisk and o for the 
text. 
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Figure 4: 
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Figure 5: 
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